Mon. Not. R. Astron. Soc. OOP. [Tlf23l f 00001 Printed 29 March 2010 (MN 1*1^ style file v2.2) 



The Dependence of Type la Supernova Luminosities on 
their Host Galaxies 

M. Sullivan^*, A. Conley^, D. A. Howell^'^, J. D. Neill^ P. Astier^, C. Balland^'^ 
S. Basa^ R. G. Carlberg^, D. Fouchez^^ J. Guy^ D. Harding I. M. Uook^'^\ 
. R. Pain^ N. Palanque-Delabrouille^^ K. M. Perrett^'^^ C. J. Pritchet^^, 
' N. Regnault^, J. Rich^^, V. Ruhlmann-Kleider^^, S. Baumont^^'^, E. Hsiao^^, 
T. KronborgS, C. Lidman^^ S. Perlmutter^^'i^ E. S. Walker^^^i 

^Department of Physics (Astrophysics), University of Oxford, DWB, Keble Road, Oxford OXl 3RH, UK 
^Department of Astrophysical and Planetary Sciences, University of Colorado, Boulder, CO 80309-0391 , USA 
, ^Las Cumbres Observatory Global Telescope Network, 6740 Cortona Dr., Suite 102, Coleta, CA 93117, USA 

• ^Department of Physics, University of California, Santa Barbara, Broida Hall, Mail Code 9530, Santa Barbara, CA 93106-9530, USA 
' ^California Institute of Technology, 1200 E. California Blvd., Pasadena, CA 91125, USA 

^ LPNHE, Universite Pierre et Marie Curie Paris 6, Universite Paris Diderot Paris 7, CNRS-IN2P3, 4 place Jussieu, 75252 Paris Cedex 05, France. 
, University Paris 11, Orsay, F-91405, France 

t ^LAM, Pole de I'Etoile Site de Chateau- Gombert, 38 rue Frederic Joliot-Curie, 13388 Marseille Cedex 13, France 

' ^Department of Astronomy and Astrophysics, University of Toronto, 50 St. George Street, Toronto ON M5S 3H4, Canada 

\ ^'^CPPM, Aix-Marseille Universite, CNRS/IN2P3, 13288 Marseille Cedex 9, France 

, INAF - Osservatorio Astronomico di Roma, via Frascati 33, OOO4O Monteporzio (RM), Italy 

• CEA/Saclay, DSM/Irfu/Spp, 91191 Gif-sur-Yvette Cedex, France 

' Network Information Operations, DRDC Ottawa, 3701 Carling Avenue, Ottawa, ON, KIA OZ4, Canada 

^'^ Department of Physics and Astronomy, University of Victoria, P.O. Box 3055 STN CSC, Victoria BC VST 1M8, Canada 
] ^^LPSC, CNRS-IN2P3, 53 rue des Martyrs, 38026 Grenoble Cedex, France 

1 Lawrence Berkeley National Laboratory, Mail Stop 50-232, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA 
' ^'^ Anglo- Australian Observatory, P.O. Box 296, Epping, NSW 1710, Australia 

' ^'^ Department of Physics, University of California, 366 LeConte Hall MC 7300, Berkeley, CA 94720-7300, USA 
, Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy 



29 March 2010 



© 0000 RAS 



2 M. Sullivan et al. 



ABSTRACT 

Precision cosmology with Type la supernovae (SNe la) makes use of the fact that 
SN la luminosities depend on their light-curve shapes and colours. Using Supernova 
Legacy Survey (SNLS) and other data, we show that there is an additional dependence 
on the global characteristics of their host galaxies: events of the same light-curve 
shape and colour are, on average, 0.08 mag (~ 4.0(t) brighter in massive host galaxies 
(presumably metal-rich) and galaxies with low specific star- formation rates (sSFR). 
These trends do not depend on any assumed cosmological model, and are independent 
of the SN light-curve width: both fast and slow-declining events show the same trends. 
SNe la in galaxies with a low sSFR also have a smaller slope ("/3") between their 
luminosities and colours with ~2.7cr significance, and a smaller scatter on SN la Hubble 
diagrams (at 95% confidence), though the significance of these effects is dependent on 
the reddest SNe. SN la colours are similar between low-mass and high-mass hosts, 
leading us to interpret their luminosity differences as an intrinsic property of the SNe 
and not of some external factor such as dust. If the host stellar mass is interpreted as 
a metallicity indicator using galaxy mass-metallicity relations, the luminosity trends 
are in qualitative agreement with theoretical predictions. We show that the average 
stellar mass, and therefore the average metallicity, of our SN la host galaxies decreases 
with redshift. The SN la luminosity differences consequently introduce a systematic 
error in cosmological analyses, comparable to the current statistical uncertainties on 
parameters such as w, the equation of state of dark energy. We show that the use 
of two SN la absolute magnitudes, one for events in high-mass (metal-rich) galaxies, 
and one for events in low-mass (metal-poor) galaxies, adequately corrects for the 
differences. Cosmological fits incorporating these terms give a significant reduction in 
(3.8-4.5cr); linear corrections based on host parameters do not perform as well. We 
conclude that all future SN la cosmological analyses should use a correction of this 
(or similar) form to control demographic shifts in the underlying galaxy population. 

Key words: supernovae: general - cosmology: observations - distance scale. 



1 INTRODUCTION 

As calibrateable standard candles. Type la supernovae (SNe 
la) provide a direct route to understanding the nature of 
the dark energy that drives the accelerated expansion of the 
Universe. Yet, the relationships that allow the calibration 
of their peak luminosities, and hence permit their cosmo- 
logical use, remain purely empirical. Relations bet ween the 
width of the SN la light curve and peak luminosity (|Phillipsl 
Il995 ) and between the SN la opti cal colours and luminosity 
fe.g. lRiess et al.lll996l :l Trippfigg^'l reduce the sca t ter in their 
peak magnitudes to ~0.1 5 mag " ijjhaet al.ll2"007l : lGuv et al.l 
I2OO7I : IConlev et aLll2008l l. As the available SN la samples 
increase in both size and quality, and the dark energy con- 
straints they provide become correspondingly more statisti- 
cally precise, it is increasingly important that the validity of 
these calibrating relationships is robustly examined. 

The observed properties of SNe la are known to cor- 
relate with the physical parameters defining their host 
galaxy stellar populations. SNe la are more than an or- 
der of magnitude more common (per unit stellar mass) in 
actively star-forming or morphologic ally late-type galaxies 
than in p assive or elliptical systems (|Mannucci et al.ll2005l : 
ISullivan et al. 2006). SNe la in elliptical or passively evolv- 
ing systems are also intrinsically fain ter, with narrower , 
faster (or lower "stretch"), light curves (|Hamuv et al.ll 19951 . 
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|l996bl: iRiess et al.lll999l : iHamuv et al]|200Gl : ISullivan et all 
200(1). Though this effect is corrected for by the light-curve 
shape correction, the amount of star formation activity in 
the universe increases with redshift, and these differences 
lead to an observed "demographic shift" in mean SN la 
properties. A greater fraction of intrinsically luminous, wider 
light-curve events in the distant universe are seen compared 
to that observed locally (|Howell et al.|[2007l ). These photo- 
metric differences are also partially reflected in SN la spec- 
tra, with SNe la in spiral galaxies showing weaker inter- 
mediate mass element line strengths than those in ellipti- 
cal galaxies (jBronder et al. lllpOS; Balland ct al. 2009), and 
a corresp onding evolution in the mean SN la spectrum with 
redshift l|SuIlivan et al.ll2009l '). 

There are suggestions that these effects may be the 
result of multiple astrophysic al channels capable of pro- 
ducing SN la explosio ns (e.g. IScannapieco fc Bildstenll2005l : 
iMannucci et al.l I2OO6I ). In particular, delay-time distribu- 
tions with distinct "prompt" and "delayed" components, 
or with a wide range of delay-ti mes, match most obser- 
vational datasets we l l (iMannucc ^ t al. 2 006: Sulliv an et al.l 
I2OO6I : IPritchet et al.ll2008l : iTotani'et all boosj ). though the 
mini mum age for the prompt systems remain s controver- 
sial (|Aubourg et al.l I2OO8I : iRaskin et all |2009| ) with some 
evidence that "prom pt" SNe la occur m ore frequently in 
metal-poor systems (jCooper et al.l [20091 ). The use of SNe 
la as precision cosmological probes therefore depends on 
establishing that the demographic shifts, or existence of 
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multiple channels to a SN la, do not impact on the light- 
curve- width/colour/luminosity relationships. If these rela- 
tions show environmental dependence, then the task of cal- 
ibrating SNe la for cosmology bec omes substa ntially more 
challenging (e.g. ISar kar et alji2008l : iKellv et al.' 2009). 

A second complication arises from the (poorly under- 
stood) colour corrections applied to SN la luminosities. Red- 
der SNe la appear fainter than their bluer counterparts, 
but the slope of the relationship between SN la colour and 
magnitude is inconsistent with the ratio of total-to-selective 
absorption appropriate for the diffuse interstellar medium 
of the Milky Way {Rv = 3.1). Multiple studies of differ- 
ent SN la samples indicate that the effective Ry inferred 
from normal SNe la is smaller than 3.1 (e.g., iTripd [l99^ : 
lAstier et al.ll200^ : iKrisciunas et al.ll2006l V and the assump- 
tion of Rv = 3.1, even after light curve shape corrections, 
leads to serious sys tematic error problems such as a sp urious 
"Hubble Bubble" (|jha et al .1120071 : 1 Conlev et al.ll2007l ). The 
reason for this low effective Rv is not well understood. Al- 
though uncorrected intrin sic variations in th e SN la popula- 
tion could play a role fe.g. lKasen et al]|2009l 'l. some dust ex- 
tinction must also affect the SN la luminosities and colours, 
and this may vary by environment. Furthermore, the exact 
value of Rv obtained is sensitive to method used to deter- 
mine it, with lower Rv obtained when fitting linear relations 
between SN la lumin osities, colours, and light curve widths 
iFolatelh et al.ll2010l ). perhaps due to intrinsic variation in 
SN la colours that correlates with luminosity but not light- 
curve width. Current knowledge of SNe la is not sufficient 
to separate and correct for both intrinsic colour-luminosity 
and dust-induced colour-luminosity effects in cosmological 
SN la samples. 

Examining how SN la luminosities vary with envi- 
ronment after light curve shape and colour corrections 
can place constraints on the degree of these possible 
variations. Early studies showed little evidence that cor- 
rected SN la l uminosities varied with host galaxy mor- 
phologies (e.g.. iPerlmutter et all Il999l: iRiess et all Il999l: 



Sullivan et alfl2003l : IWilliams et al.ri2003l : [Gallagher et all 
20051 ). though these tests used relatively small samples of 



events (< 50), in some cases from the first-generation of 
SN la cosmological samples before dense multi-colour light 
curves were routinely obtained. 

More recent analyses, using larger, well-observed 
samples, have show n tentative evidence for variation. 
iHicken et al] (|2009bD found ~ 2a evidence that SNe la in 
morphologically E/SO gala^xies are brighter than those in 
later-type spirals after light-curve shape and colour correc- 
tions. Extending beyond simple host galaxy morphologies to 
more physically motivate d variables gives f urthe r tantalising 
suggestions of variation. iGallagher et al.l (|2008l ) found evi- 
dence for a correlation between Hubble diagram residual and 
host galaxy stellar metallicity in a sample of 17 local SNe 
la located in E/SO galaxies, in the sense that fainter SNe 
la after correction were found in metal poor systems (note 
this is the reverse of the originally published trend due to 
an error in t he original analysis; P. Garnavich, private com- 
munication). iHowell et al. I (20091 ) used 55 SNe la from the 
first year of the SNLS and showed no significant correla- 
tion between Hubble residual and host galaxy metallicity, 
albeit using host gas-phase metallicities inferred from av- 
erage galaxy stellar-mass-metallicity relations, a less direct 



measure of metallicity. iKellv et all (|2009l ) have shown a re- 
lation between host galaxy stellar mass and Hubble residual, 
in the sense that more massive systems host brighter SNe 
la after light curve shape and colour corrections. Under the 
assumption that more massive galaxies are metal rich, thi s 
trend is consistent with the revised Ga llagher et al.l (|2008l ) 
result. 

In this paper, we use a sample of 282 high redshift SNe 
la discovered and photometrically monitored by the Canada- 
France-Hawaii Telescope (CFHT) as part of the Supernova 
Legacy Survey (SNLS), and which form the SNLS "three- 
year" sample (SNLS3). Using deep optical imaging of their 
host galaxies taken over the duration of the survey, we place 
constraints on their recent star-formation activity, stellar 
masses (and hence inferred metallicity), and compare to the 
photometric properties of the SNe la that they host. In par- 
ticular, we search for evidence that the corrected SN la lumi- 
nosities correlate with these host properties, indicating pos- 
sible systematic errors in the light curve fitting framework 
that underpins their cosmological use. We compare with the 
properties of a sample of lower-redshift SNe la taken from 
the literature. 

A plan of the paper follows. In §[2] we introduce the SN 
la sample and the data available on their host galaxies. § O 
investigates how the SN la light curve widths and colours of 
these SNe la varies according to their host galaxy proper- 
ties, and in §|3]we compare their corrected luminosities to the 
host properties. We discuss the results, including the cosmo- 
logical implications, in §0 and conclude in §[H1 Throughout, 
where relevant we assume a fiat ACDM cosmological model 
with JIm = 0.256 (the reason for this non-standard choice 
is explained in § [3)) and Ho=70kms^^ Mpc^^ assumed in 
all quo ted abso l ute m agnitudes. All magnitudes are given 
on the Landoltl ll 19921) photometric system as described in 
iRegnault et all l|2009l ). 



2 TYPE lA SUPERNOVA AND HOST 
GALAXY DATA 

We begin by introducing the SN la samples that we will use 
in this paper, and the associated data available for their host 
galaxies. 



2.1 The SN la samples 

The high-redshift SN la data is taken from the Super- 
nova Legacy Survey (SNLS). This used optical imaging data 
taken as part of the deep component of the five-year Canada- 
France-Hawaii Telescope Legacy Survey ( CFHT-LS) us- 
ing t he square-degree "MegaCam" camera l|Boulade et al.l 
|2003| ). located in the prime focus environment "MegaPrime" 
on the CFHT. The "deep" component of CFHT-LS con- 
ducted repeat imaging of 4 low Galactic-extinction fields, 
time-sequenced with observations conducted every 3-4 
nights in dark and grey time. Four filters, gMrMtMZM, 
were used allowing the construction of high-quality multi- 
colour SN light curves; um data were also taken but are 
not time-sequenced. On each night of observations, the data 
were searched using two independent pipelines, and an amal- 
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gamated candidate list producecQ fsee lPerrett et aLll201Cf ). 
Spectroscopic follow-up confirmed SN types and measured 
redshifts to allow SNe la to be placed on a Hubble diagram. 

In this paper we use SNe la belonging to the three- 
year sample, SNLS3; this includes all SNe la discovered up 
until the end of July 2006. The lig ht curves and oth er de- 
tails for these SNe can be found in IGuv et all (l2010l'l. ari d 
spectroscopic i nf orma tio n is taken fr o mllfowell et al.l II2005I). 
iBronder et all (|2008l ). lEllis et all (j2008h . iBalland et all 
I 2009( ). and Walker et al. (2010). The fuU SNLS3 sample 
consists of 282 SNe la, however we exclude some of these 
events from some of our analyses: 



(i) Only SNLS SNe la with an identified (§ [O]) host 
galaxy (272 events) are consid ered (for a discussion of the 
identification procedures, sec iSuUivan et al, .2006 , hereafter 
S06), 

(ii) Only SNe la passing light c urve quality cuts (e.g. 
IConlev et aTll2008l : IGuv et al.|[2010l ') are used - there must 
be sufficient photometric coverage to derive reliable peak 
luminosities, lig ht curve widths, and colours (see details in 
IGuv et al.|[2oTol ). 

(iii) We only consider normal SNe la with light-curve pa- 
rameters in the range considered for cosmological fits - our 
motivation in this paper is to assess the cosmological im- 
pact of any host galaxy dependent trends. Specifically, we 
require that the stretch of the SN be greater than 0.80, and 
that the colour be less than 0.30 (see S 12.21 for a discussion 
of the meaning of these parameters). We also discard > 3(t 
outliers from the best fitting cosmological model. 231 events 
pass the light-curve coverage and SN parameter cuts. 

(iv) In analyses where we search for trends in the SNLS 
data, we use only SNe la at redshift z ^ 0.85 (195 events). 
At these lower redshifts, both the SN and host galaxy pho- 
tometry are higher signal-to-noise and their photometric pa- 
rameters bet ter measured. The SNLS Malmquist biases are 
also smaller (jPerrett et al.ll201Gl l. 

Where relevant, we also use samples of SNe la from 
the literature. We const ruct a sample of low -redshift SNe la 
from the compilation of lConlev et al.l l|2010l) , which includes 
SNe la from a variety of sources (primarilv | Hamuy et al.l 



1996a'; 'Ricss ct al." 19991; |jha et al.ll2006l : iHicken et al.l(2009al : 

Contrcras ct al. 2009). We apply bulk-fiow peculiar veloc- 
ity corrections to the SN magnitudes and redshi fts, placing 
the re dshifts in the CMB-frame (z cmb) following Neill et al] 
|2003), but with updated models (|Conlev et al.ll2010l '). The 
accuracy of these corrections is estimated to ±150 km s~^, 
which we propagate into the SN magnitude errors in cosmo- 
logical fits. We only use SNe la in the smooth Hubble flow, 
here defined as Zcmb ^ 0.01, and apply the same light curve 
quality cuts as for the SNLS sample. There are 110 low- 
redshift o bjects in total. W e also use the ffST-dis covered 
sample of iRiess et al.l (|2004l ') and iRiess et al.l (|2007l '). here- 
after the R07 sample. We select 24 SNe la at z > 0.9 from 
this sample to increase our redshift lever arm above z = 1. 



^ Candidates can be 

http: / /legacy . astro .utoronto . ca/ 



found 



at 



2.2 SN Light curve fitting 

In its current application, SN la cosmology depends on two 
corrections to "raw" SN la peak luminosities that when 
applied reduce the dispersion in their peak magnitudes. 
The first is t he light-curve-shape/luminosity relation (e.g. 
|Phillips|[l99^ ): brighter SNe la tend to have wider, longer- 
duration light curves (higher stretch) than their fainter 
counterparts. The second is a colour correction: brighter SNe 
la tend to have bluer colours, whilst fainter SNe tend to be 
redder l|Riess et al.ll 19961 : iTripdl 19981 ). Together the applica- 
tion of these corrections can yield distance estimates precise 
to ~ 6 per cent. These corrections are applied in different 
ways depending on whet her a technique is a distance es- 
timator (e.g., MLCS2k2; |jha et al.l 120071') or a light curve 
fitter (e.g. SALT or SiFTO: IGuv et al.ll2007l : IConlev et all 
'2008)), though the underlying principle is the same in both 
approaches. 

I n this paper, we pr imarily use the SiFTO light curve 
fitter (IConlev et al.ll2008l ') and compare our results to SALT2 
(IGuv et al.l 2007 ) where appropriate. In general SALT2 and 
SiFTO give very similar results when trained o n the same 
SN sa mple - a full discussion can be found in IGuv et al] 
(|2010l ). Both fitters have been retrained and improved since 
the original published versions using SNLS and low-redshift 
data. The product of both fitters for each SN is a rest-frame 
_B-band apparent magnitude (m^), a stretch (s) measure- 
ment, and a colour estimate (C), together with associated 
errors and covariances (SALT2 uses the broadly equivalent 
x\ parameter in place of s). Throughout, the SN la colour, 
C, is defined as the rest-frame {B — V) colour of the SN at 
the time of maximum light in the rest-frame B-band. We 
refit all SNe la using these light curve fitters to ensure that 
the different samples can be placed on the same distance 
scale. A full discussion of their application to the current 
d ataset, together w ith t heir tabulated outpu t, can be found 
in lGuv et al.l (|2010l ) and lConlev et all (|2010l '). 



2.3 SN liost galaxy photometry 

Our SNLS host galaxy photometry comes in the optical from 
the CFHT-LS {uMQAirMiMZM) and in the near infra red 
(IR) from the WIRcam Deep Survey (WIRDS; Bielby et al. 
in prep.) of a sub-section of the CFHT-LS fields (J, H, Ks). 
The identification proced ure fo r the SNLS SN la host galax- 
ies is the sa me as that in | S06l. Photom etry is performed by 
SExtractor (|Bertin fc ArnoutsI [l996l ) using FLUX_AUTO 
photometry running in dual-image mode, detecting from the 
deep iM stack and measuring from each of the optical and 
near-IR filters in turn (60% of our SN la hosts have near-IR 
data). Each stack has a similar image quality and hence the 
same physical aperture is used in each filter. In about ~3% of 
cases no host galaxy can be identified. This could be because 
the SN lies far from the centre of its host galaxy leading to 
ambiguity in the correct choice of host, or because the host 
lies below the CFHT-LS fiux limits. We discard these ob- 
jects. Weight maps are used to determine the measurement 
errors, and in the optical, the photometric zeropoints are 
generated usin g a compariso n to the tertiary standard star 
lists of Rogna ult et al.1 ()2009h . No SN light is present in the 
optic al stacks which are constructed on a per season basis 

ilsoeh. 
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For the low-redshift s ample, we u s e host galaxy photom- 
etry recently compiled bv lNeill et~all (|2009l . hereafter N09), 
including ultraviolet data from the GALEX (Galaxy Evo- 
lution Explorer) satellite, and optical photom etry from the 
third reference catalog of bright galaxies fRC3: [Corwin et alj 
119941 ) and the Sloan Digital Sky Survey (SDSS). The pho- 
tometry for this sample was carefully performed in matched 
apertures with foreground contaminating stars masked. 
Though these data span a different observed wavelength 
range compared to the SNLS host photometry, in the rest- 
frame the wavelength range covered is reasonably similar: 
the maximal range is 1500-9000 A for the low-redshift sam- 
ple (though frequently the available data span a smaller 
range than this), and 2400-13000 A for the SNLS sample (at 
z — 0.6). Note that not all the low-redshift SN hosts have 
publicly available host photometry: only 69 (of 110; 63%) 
low-redshift SNe la have sufficient data and are carried for- 
ward in the analysis. The missing low-redshift SNe la are 
due to incomplete GALEX and SDSS coverage, rather than 
the host galaxies being too faint to be detected by the two 
surveys. 

For the R07 sample, we use photometry taken from 
the Great Observatories Origins Deep Survey (GOODS) 
HST data (Gia valisco et al. 2004), taken with the Advanced 
Camera for Surveys (ACS). Data are available in four fil- 
ters: F435W (broadly equivalent to a B filter), F606W (V), 
F814W (i') and F850LP {z'). All of the R07 SNe have ACS 
coverage, and again SExtractor FLUX_AUTO photometry 
is used. We also use J, H and K imaging taken as par t of 
the GOODS NICMOS surv ey (e.g.lBuitrap et alll2008l ). as 
well as ground-based data 



r ev le.g. I -Duitraeo et a 
(|Retzlaff et al.ll2009l V 



2.4 Host galaxy parameter estimation 



of 



The method to estimate physical parameters 
the SN la host galaxies is similar to that in 
SOSi wh ich used the photometric re d shift code Z- 
PEG l|Le Borgne fc Rocca-Volmerangd I2OO2I ) based 
upon the PEGASE.2 s pectra l synthesis code (e 



iFioc &: Rocca-Volmerang^ 1 19971 ). We use an expanded 
set of 15 exponentially declining star formation histories 
(SFHs) with SFR oc exp"*''^, where t is time and r is 
the e-folding time, each with 125 age steps. The internal 
PEGASE.2 dust prescriptions are not used, and instead a 
foreground dust screen varying from E{B — V) = to 0.30 
in steps of 0.05 is added. With the 7 different foreground 
dust screens, this gives a total of 105 unique evolving 
spectral energy distributions (SEDs). The metallicity of 
the stellar population evolves consistently following the 
standard PEGASE.2 model with an initial value of 0.0004, 
and the standard PEGASE.2 nebular emission prescription 
is added. 

Z-PEG is used to locate the best-fitting SED model 
(in a sense), with the redshift fixed at the CMB- frame 
redshift of the SN. Only solutions younger than the age of 
the Universe at each SN redshift are permitted. The current 
stellar mass in stars (Mgteiiar, measured in Mq), the recent 
star formation rate (SFR, in MQyr~^, averaged over the last 
0.25Gyr before the best fitting time step), and the specific 
star formation rate (s SFR, the SFR per unit stellar mass 
with units of yr~^, e.g. lCuzman et al]|l997l ) are all recorded. 
Error bars on these parameters are taken from their range in 



the set of solutions that have a similar (^s in ,S06i ). Note 
that we do not measure the instantaneous SFR as we only 
fit broad-band photometry. We refit all the [nO^ photometry 
to ensure the e xact same libra r y SED s are used for all hosts. 

We use a i Ran a fc Bas ul (Wdf) initial mass function 
(IMF), the PEGASE.2 default, throughout. Our results are 
not sensitive to this choice - we have repeated our a nalysis 
in full wit h both the more standard Salpeteij (|l955l ) IMF, 
and with a iBaldrv fc Glazebrookl (|200a 7 IMF, and find sim- 
ilar results, though the Mgtoiiar and SFRs derived for the 
host galaxies have (expected) small mean offsets when using 
different IMFs. In detail, the use of a Salpeter IMF gives sys- 
tematically larger host masses by 0.04 dex (smaller masses 
by 0.16 dex for the B&G IMF), and smaller SFRs by 0.04 
dex (smaller by 0.16 dex for B&G), with scatter of around 
O.ldex in each comparison. These differences are not a func- 
tion of A/stciiar, SFR or redshift and so have a negligible 
impact on our conclusions. 

As only ~60% of our SNLS SN la hosts have observer- 
frame near-IR data, we compare the derived properties with 
and without these data to check for potential biases in the 
remaining 40% of objects. The mean difference in Mgteiiar 
(defined as M°^^i^^-M™iiar) is 0.001 dex (r.m.s. 0.15) and 
in the recent SFR the mean difference (SFR°'"^-SFR™) is 
—0.18 dex (r.m.s. 0.44), in the sense that excluding the IR 
data leads to smaller SFRs. Thus we find no evidence that 
the near-IR data leads to systematically different Mstoiiar, 
and mild evidence that including these data leads to larger 
SFRs. The differences in SFR do not follow a Gaussian dis- 
tribution; instead the difference is centred around zero but 
with a long tail to negative differences; therefore we choose 
not to apply the mean offset to the 40% of hosts with no 
IR data. There is no evidence for any redshift dependent 
trend. Information on the derived properties for the SNLS, 
low-redshift and R07 hosts can be found in Table [T] 

The MstcUar and SFR distributions for the SNLS and 
low-redshift samples are shown in Fig. [1] together with the 
distribution of galaxies in the SFR-Msteiiar plane. As might 
be expected, galaxies with the smallest sSFRs tend to be the 
most massive systems, with the lowest mass systems univer- 
sally consistent with strong star formation activity. As pre- 
viously highlighted bv lN09l . the SNLS and low-redshift hosts 
show quite different distributions in Mstdiar (and to a lesser 
extent SFR): The low-redshift SNe are drawn from more 
massive host galaxies. This is almost certainly due to selec- 
tion biases. SNLS is a rolling search which will locate SNe 
la in any type of host galaxy in which they explode, and, 
modulo any small spectroscopic follow-up bias, this range 
will be reflected in the cosmological sample. At low-redshift, 
most SNe la are drawn from galaxy-targeted searches which 
search known (and typically bright/massive) galaxies, conse- 
quently the m ost massive systems will be over-represented. 

Following iHowell et al] l|2009l ). we convert the Afstciiar 
mass estimates into metallicities using average Mstdiar- 
metallicity (Mstoiiar-^) relations. As the universe ages, 
galaxies will become more massive via merging processes, 
and more metal rich following chemical enrichment and de- 
creased metal loss. We use a relation between gas-phase 
metallicity, explicitly the nebular oxygen abundance relative 
to hydrogen, O/H, and A fstciiar derived from SDSS galax- 
ies l|Tremonti et al.ll2004h . We use units of 12 + log(0/H), 
where the solar abundance is ~8.69. This MstcUnv-Z rela- 
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Table 1. Properties for the SNLS, low-redshift and ROT host galaxies. The full table can be found in the electronic version of the journal. 



SN Name ^cmb *Af LOG A4tciiar LOG SFR 

(Mq) (Mgyr-l) 
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-0.77 
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q 55+0.13 
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+0.58 
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Figure 1. The distribution of SN la host galaxies in Mstollan SFR and inferred gas-phase metallicity. Left: Histograms in Mstollan SFR, 
and metallicity; low-redshift (open histogram) versus SNLS (filled histogram). Right: The distribution in the Mateliar~SFR plane. The 
dashed line is a line of constant sSFR used as the default split to separate high sSFR and low sSFR galaxies. The light-grey shaded areas 
show the range over which the default splitting value is varied in Section 14.21 SNLS SN la hosts are shown as filled circles, low-redshift 
hosts as white circles, R07 hosts as grey circles. Note that the apparent diagonal ridgeline in the Mg^^nj^j-SFR plane is a result of the 
maximum sSFR being reached in the model SFHs, a disadvantage of the simple SFHs considered here. 



tion is observed to evolve with redshift, being shifted to- 
wards high Afst ciiar or lower metalliciti es at higher red- 
shift (e.g. Savaglio et aLllioOsI : [Lamarcil le et al.ll2009l 'l. To 
acc ount for this ef f ect, w e use the evol ution measured 
by iLamareille et all (|2009l ') relative to the iTremonti et al.l 
relation. (Note that the use of a Msteiiar--^ re- 
lation that evolves with redshift will exacerbate the dif- 
ference in Matciiar between low redshift and SNLS hosts 
when expressed as a metallicity.) Though the calibration of 
nebular-line metallicity diagnostics is a complex topic (e.g. 
iKewlev k, Ellisonll2008l ). the exact calibration scale does not 
concern us here, so long as we can place all our hosts on the 
same relative system to search for variations in the SN sam- 
ple. 

The ITremonti et al.l l|2004 ) relation is derived for gas- 
phase metallicity, and may not be applicable in elliptical 
galaxies with little cold gas where a stellar metallicity may 
be more appropriate. In principl e, we could also use a stellar 
metallicity-mass relation (e.g.. [Callazzi et al.ll2005l ). Ifow- 
ever, this relation has only been measured at low-redshift, 
and any (expected) evolution with redshift has not been 



constrained observationally - stellar metallicities are signif- 
icantly more difficult to measure than gas-phase metallici- 
ties, requiring absorption line measures instead of emission 
lines. Therefore we restrict ourselves to the gas-phase metal- 
licity, but note that as the mechanism for retaining metals 
is the depth of the galaxy gravitational potential well, el- 
lipticals (or passively evolving galaxies) should follow the 
same general trend between stellar metallicity and A/gteiiar 
(|Gallazzi et al.|[2005l ). with stellar met allicity and gas-phase 
nebu l ar metallicity well co rrelated (e.g. lCid Fernandes et all 
l2005l : iGallazzi et aLlbOOSl ) . 

The broad-band SED fitting approach used here is a rel- 
atively crude way to determine galaxy properties. However, 
our choice is limited as, for the SNLS sample at least, we 
only have optical UMginrMiMZM magnitudes (albeit mea- 
sured from extremely deep and well-calibrated images) and 
some near-IR information for the hosts. In particular, spec- 
troscopic diagnostics used in many other galaxy analyses are 
not available. Though each SN la does have a spectrum con- 
taining some host galaxy spectral information, these form a 
very heterogeneous set invariably contaminated by SN light. 
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making the accurate measurement of host spectral features 
impossible. We did explor e add ing GALEX ultraviolet data 
to our SNLS host fits (as inlN09ll: however the relatively poor 
resolution of the GALEX imaging led to substantial source 
confusion for most of our galaxies. We also experimented 
with the use of more complex ( "stochastic" ) SFHs involving 
bursts of star formation superimposed on smooth underlying 
SFHs as in ot her studies of SN la or gamma-ray bur st host 
galaxies (e.g. ISavaglio et ai]|2009l : ISchawinskil |2009| ); how- 
ever the increased number of free parameters (e.g., burst 
strength, burst duration, burst age, etc.) makes the prob- 
lem quite degenerate. Hence, our host galaxy parameters are 
only representative of the dominant star formation episode 
in each host. We note that the use of stochastic SFHs re- 
sults in i ncreased derived stellar masses of galaxies by about 
0.14 dex jPozzetti et al.ll2007l : lLamareille et ahlbood v 

Where required in our analysis to examine the depen- 
dence of SN la properties on environmental properties, we 
classify the SN la host galaxies into different groups. The 
first split is performed according to their sSFR: galaxies with 
log(sSFR) < —9.7 with smaller amounts of star formation 
relative to their stellar masses are classified as low-sSFR, 
and those with a larger sSFR are classified as high-sSFR. 
This split will separate those galaxies dominated by young 
stellar populations from those comprised of older stellar pop- 
ulations. The second split is based upon the host galaxy 
A/stellar (log(Mateiiar ) > 10.0 are defined as high mass), with 
an extension to the inferred gas-phase metallicity {Z > 8.8 
are defined as metal rich). (We consider the effect of varying 
these split points in later sections.) Approximately 30 per 
cent of SNe la in our SNLS sample are found in low-sSFR 
systems according to the definition above. 



3 SN lA PROPERTIES AS A FUNCTION OF 
HOST GALAXY PROPERTIES 

In this paper we correct SN la peak magnitudes for the 
stretch and colour relations (m§"^'^) using a standard empir- 
ical approach with the form: 



mg'^'^ = rriB + a [s — 1) — pC 



(1) 



where a parametrizes the stretch-luminosity relation, /3 
parametrizes the colour-luminosity relation, and tub, s and 
C are the observed peak magnitude, stretch and colour out- 
put from the light curve fitters (i; 12. 2|) . ms and C are cor- 
rected for Milky Way extinction but not for any host galaxy 
extinction, a and /3 are typically determined fro m simultane- 
ous fi ts with the cosmological parameters (e.g. lAstier et al.l 
I2OO6I . see also §2]), and reduce the scatter in the peak SN la 
magnitudes to ~0.15mag. 



3.1 Light curve shape and colour 

The correlation between SN la light curve s hape (stretch, 
Ami s , etc.) and host galaxy mo r phology (e.g. Hamuv et al.l 
I199E] . Il996bl : iRiess et all 1 19991: iHamuv et all [2OO0I ). SFR 
l S0g : lN0g| l, and Mateiiar l|Howell et all 120091 ') is well estab- 
lished: fainter, lower stretch (high Amis) SNe la explode in 
older stellar populations showing little ongoing star forma- 
tion. 



We extend this earlier work to examine the SN la prop- 
erties as a more continuous function of sSFR and Mstoiiar in 
Fig. [21 The expected trends with stretch are recovered. Low 
stretch SNe la preferentially explode in low-sSFR galaxies 
with little or no ongoing star formation, and are rarer in 
galaxies which have substantial star formation. By contrast 
the higher stretch SNe la are found in galaxies with a range 
of mean sSFRs but with a deficit in the low sSFR systems. 
Similar trends are seen as a function of Msteiiar: lower-stretch 
SNe are almost entirely absent in low-Mstdiar galaxies (ex- 
pected to have the highest sSFRs). There is is also evidence 
for stretch being a continuous variable of sSFR or Mstciiar- 

No conclusive trends between SN la rest-frame colour 
and host galaxy properties have yet been identified. In our 
data, there are also no strong trends, though SNLS SNe la 
in low sSFR systems do show slightly bluer colours in the 
mean (no differences are seen between the colours of SNe la 
in low-Afsteiiar hosts and high- Mstoiiar hosts). A t-test shows 
the mean colours of SNe in low and high sSFR galaxies have 
about a 97% chance of being different; a Mann- Witney U 
test gives a similar 97% chance that the colours arise from 
different distributions. This is consistent with the simple 
viewpoint that more strongly star-forming galaxies contain 
more dust, and that dust will make SN colours redder. How- 
ever, the lack of a large difference in mean colour suggests 
that the amount of dust along the line of sight to SNe la in 
star-forming host galaxies is small. 

Note that the lack of a trend in the opposite sense may 
also be considered surprising at first sight. There is a strong 
variation in stretch with star formation activity, and ob- 
servations show that the very lowest-s (and typically sub- 
luminous) SNe la are usually redder in C - this would imply 
that the lowest sSFR systems should host the reddest SNe, 
which is opposite to the weak trend that we do observe. 
However, the strength of this stretch-colour relati onship for 
"normal" SNe la at maximum lig ht is weak (e.g. |jha et al.l 
l2006l : lNobih fc Goob ar 2008; Tak anashi et al.ll2008l 'l. Trends 
do exist between oth er colours and stretch and at other light- 
curve phases (e.g. iPhillips et al. I I1999I : iNobili fc Goobaj 
,200 i) and for sub-luminous SN1991bg-li ke very low-s SNe 
la ( Garnavich et al. |[200i: IJhaet ZI lBo^). but these are ex- 
cluded from our analysis by the s > 0.8 requirement. 



4 SN LUMINOSITY DEPENDENCE ON HOST 
CHARACTERISTICS 

We now examine the effect of the environment on the "cor- 
rected" brightness of SNe la, and hence the effect on their 
use in a cosmological analysis. We minimize 



(2) 



N 

where m'^g" is given by eqn. ((T)), CTstat is the total identified 
statistical error and includes uncertainties in both tub and 
m^°'^, (Tint parametrizes the intrinsic dispersion of each SN, 
and the sum is over the N SNe la entering the fit. rn^°'^ is 
the model S-band magnitudes for each SN, given by 



mB° = blog^QVL {z;Qm) + Mb, 



(3) 



where 2?^ is the c/Hq reduced luminosity distance with the 
c/Ho factor (here c is the speed of light) absorbed into Mb, 



© 0000 RAS, MNRAS 000.[Tlf23l 



8 M. Sullivan et al. 



1.3F 



£ 1.0 r 

55 



0.8 T 
0.7 L 




0.4- 



0.2 ■ 



O 

^ 0.0 ■ 



-0.2 ■ 




-12 



-11 -10 -9 

LOG host galaxy sSFR (yr"') 



-12 -11 -10 -9 

LOG host galaxy sSFR (yr ') 




9 10 11 

LOG host galaxy M3„„3,(Ms) 



9 10 
LOG host galaxy M,,,,,,, (Me 



Figure 2. The SN la stretch (s; left panels) and colour (C; right panels) from SiFTO as a function of the host galaxy sSFR (upper 
panels) and M^^^n^^ (lower panels). The red points show the weighted mean s or C, co rrecte d for dispersion, in bins of sSFR and M^^^^n^-^. 
Only SNLS SNe are shown (similar plots for the low-redshift sample can be found in lN09t) . Equivalent results are found with SALT2. 



the absolute luminosity ofas = lC = SN la (eqn. (HJ). 
Explicitly, Mb = Ms -|- 5 logi(,(c/ii"o) + 25, where Mb is the 
absolute magnitude of a SN la in the B-band (for SALT2 
fits, Mb refers to an a;i = 0, C = event). For convenience, 
we present our results as Mb rather than ^A b , but note that 
this requires a value of Hq to be assumed - a choice that 
does not impact our results in any way. 

a, P and Mb are often referred to as "nuisance vari- 
ables" in cosmological fits as they are not of immediate in- 
terest when determining cosmological models. Instead they 
parametrize luminosity variations within the SN la samples 
and are likely related to the physics of the SN la explosion 
and /or the SN la environment; clearly this makes them of 
great interest in this paper. 

Two different approaches are used. The first approach 
examines the residuals of the SNe from global cosmologi- 
cal fits using the SNLS3 z < 0.85 plus low-redshift sample, 
fixing Q,M = 0.256 (the best-fit value for this sample). We 
choose this number instead of a more "standard" value like 
Qm = 0.3 to ensure that no redshift bias in our SN la resid- 
uals is introduced by adopting a cosmological model that 
does not fit the data adequately. The second examines any 



variation of the nuisance variables by fixing the same cosmo- 
logical model and performing fits on sub-samples of SNLS 
SNe with the nuisance parameters free. The first technique 
uses global values of the nuisance variables for the entire 
sample, whereas the second allows them to vary by environ- 
ment. The advantage of the latter technique is that as the 
cosmological model is fixed, a large low redshift sample is 
not required in order to examine brightness-dependent tests 
internally within the well-measured SNLS sample. Through- 
out, we define a Hubble residual as mg"^"^— m^"'^, i.e. brighter 
SNe have negative LIubble residuals. 



4.1 Residuals from global cosmological fits 

We consider the residuals from the best fitting cosmological 
model as a function of three host properties: sSFR, Mstoiiar, 
and Mgtciiar converted into a metallicity estimate (Z). Resid- 
ual trends here indicate luminosity-dependent effects that 
are not removed by the standard s (or xi) and C methodol- 
ogy, but that do correlate with some other physical variable 
associated with the host galaxy. We emphasise that Mgteiiar 
and Z are therefore highly correlated, and our Z estimates 
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are depende nt on the evolution in the MsteUa,T-Z relation 
measured bv lLamareille et all l|2009l) . 



4.1.1 Host specific SFR 

The residuals in the SNLS and low-redshift samples as a 
function of host galaxy sSFR are shown Fig. (the R07 
sample has too few SNe for any meaningful trends to be 
detected). In the SNLS sample, there is evidence that SN la 
events in low sSFR galaxies appear brighter on average than 
those in high sSFR galaxies after s/x-i and C corrections. The 
numerical values for the mean residuals in each bin are given 
in TableElfor SiFTO and Table[3]for SALT2; the differences 
seen are not dependent on the light curve fitter used. Note 
that this trend is opposite to that perhaps expected from 
the observation that low-stretch (intrinsically fainter) SNe la 
reside in lower-sSFR hosts (S I3.1I) : however this observation 
is pre-stretch correction - our residuals have been corrected 
for both stretch and colour (this is discussed in more detail 
in §[531). 

The differences in the SNLS samples between the two 
lowest and two highest sSFR bins are ~2.6cr significance. In 
the low-redshift sample analysed here, the trends are con- 
sistent with the SNLS sample (discarding the most strongly 
star-forming bin with only a small number of events), but 
low significance. For the SNLS data, fitting a straight line 
to the binned points and accounting for errors in both axes 
detects a non-zero gradient at ~2.5cr. 



4.1.2 Host stellar mass 

The residuals as a function of host galaxy Mstoiiar are shown 
in Fig. U In this case, SNe in more massive galaxies appear 
brighter on average than those in lower-mass galaxies, af- 
ter s/xi and C corrections. The mean residuals are given 
in Tables [2] and |3l middle panel. In the SNLS sample, SNe 
in the three lowest Mstoiiar bins are fainter than those in 
the two highest Matciiar bins at ~3.9cr significance. The low- 
redshift sample is consistent with this, though again the sig- 
nificance is smaller and the range in Mstoiiar is much more 
restricted. For the SNLS data, fitting a straight line to the 
binned points detects a non-zero gradient at ~3.1ct. 



4.1.3 Host inferred metalUcity 

We convert our Mstoiiar into gas-phase metallicities [Z) as 
described in i) l2.4l The residuals as a function of host galaxy 
Z are shown in Fig. [S] As the Z estimates are directly re- 
lated to the Mstoiiar, the same trends are apparent - indeed, 
had we used an M^tciui-Z relation that did not evolve with 
redshift (e.g. Trcmonti ct al. 2004; Gallazzi o_t al. 20051). the 
Z and Mstoiiar results would be identical. The mean resid- 
uals are given in Table [5] and O lower panel. In the SNLS 
sample, SNe in \ow-Z galaxies are fainter than those in high- 
Z galaxies at ~3.7(t significance. The low-redshift sample 
is consistent with this, though the range in Z is more re- 
stricted. Fitting a straight line to the binned SNLS points 
detects a non-zero gradient at ~3.0o". 



4.2 Universality of the nuisance variables 

We now consider the s/xi and C corrections that are applied 
to the SN magnitudes, and test whether these two correc- 
tions are consistent for sub-samples of SNe divided accord- 
ing to the characteristics of their host galaxies. We perform 
these tests using the larger and more complete SNLS SN la 
sample, and again hold the cosmological model fixed. 

As discussed in Section 12.41 we separate the hosts ac- 
cording to their sSFR, Mstoiiar, and Z. The exact values 
chosen as the split points are a somewhat subjective choice 
(e.g. Fig. [TJ. The Mstoiiar split point was chosen to sepa- 
rate the hosts into bins of approximately equal sizes, and 
the Z split point is this Mstoiiar value converted into a Z at 
z = 0.5; no fine-tuning was attempted. We explore how our 
results change as we modify the default splits of § 12.41 over 
a range of 0.8 in log(sSFR), 1.0 dex in log(Msteiiar), and 0.2 
dex in Z. We perform fits for a, (3 and Mb, in each case 
iteratively adjusting aint until a per degree of freedom of 
1.0 is obtained. The SiFTO results are given in Tables g] [5] 
and|Bl 

As expected from §133] there are several host dependent 
trends. In particular, we find that 

• SNe la in high- Mstoiiar or high-Z hosts prefer brighter 
Mb than those in low-Afstoiiar (low-Z) hosts (~ 4. Oct), as do 
SNe in low-sSFR versus high-sSFR hosts (~ 2.5a), both in 
the same sense as the trends in Fig. [31 

• SNe la in low-sSFR hosts show smaller values of /3 than 
those in high-sSFR hosts (~ 2.7a). Smaller (< 2.5cr) differ- 
ences are also seen between /3 measured in high and low- 
Msteiiar (Z) hosts, though this is more dependent on the 
split point used, 

• The a values in the split samples are consistent, differ- 
ing only at < 1.5a, 

• SNe la in low-sSFR hosts show a smaller r.m.s. scatter 
about the best fits than those in high sSFR hosts; the r.m.s. 
scatter of SNe in high-A/stoiiar {Z) and low-i\fstoiiar (Z) hosts 
are similar. 

These results are not sensitive to the split points used to 
characterise low-sSFR versus high-sSFR, or high-Afstoiiar 
(Z) versus low-A/stoiiar (Z) hosts. Fig.|6]shows the joint con- 
fidence contours in the three combinations of the nuisance 
variables for the sSFR and Mstoiiar split samples. 

We also assess the significance of these results using 
a Monte-Carlo permutation test of our data. We randomly 
draw A^iow-sSFR (or iVhigh-M when considering the Afstoiiar 
split) SNe la from the full SNLS sample (without replace- 
ment) to generate a fake low-sSFR (high-Afstoiiar) sample 
with the same number of events as the true low sSFR sam- 
ple. The A'high-sSFR (A^Iow-m) remaining SNe are assigned 
to the fake high sSFR (low-Mateiiar) sample. The nuisance 
variable fits are repeated on both datasets, and the best- 
fitting Q, P and Mb recorded. We repeat this procedure 
25000 times, and compare the distribution of derived nui- 
sance variables from the permuted samples to the values 
recovered for the real samples. We assess how frequently the 
difference between the nuisance variables for the fake sample 
is equal to or greater than the difference between the same 
parameters measured on the real data. The Monte Carlo 
approach can also assess the significance of any differences 
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Figure 3. Residuals from the best-fitting cosmology for the SNLS (top) and low redshift (bottom) samples as a function of host galaxy 
sSFR. Hosts with s < 1 SNe la are shown as open circles, and those with s ^ 1 as filled circles. Brighter SNe la after correction have 
negative residuals. The error bars on the individual SNe are taken from the SED fitting for the sSFR axis, and are the statistical errors 
propagated through the light curve fitting for the residual axis (i.e. they do not include the crint component). The red circles are the 
mean residuals in bins of host sSFR, drawn at the mean value of the sSFR in each bin. The error bars on these points represent the bin 
width in sSFR, and the weighted error on the mean (corrected for dispersion in each bin) for the residual axis. The right-hand histograms 
show the total residuals for SNe la in low and high sSFR hosts. 



seen in r.m.s. scatter between different samples, otherwise 
difficult to assess statistically. 

For samples split by sSFR, we find that the differences 
in a are found in randomly permuted samples 23% (a ~ 
1.2a significance) of the time, for /3 2.6% (~ 2.2o-), for Mb 
0.5% (~ 2.8a) and for the r.m.s. scatter 3.9% (~ 2.1a). For 
samples split by Msteiiar , the differences are found in random 
samples 34% (~ 0.95a) for a, 6.5% (~ l.Scr) for /3, 0.005% 
(~ 3.5cr) for Mb, and 44% (~ O.Scr) for the r.m.s.. These 
results support our main finding above: the most significant 
differences are seen in Mb for all methods of splitting the 
host galaxies. Monte Carlo tests show smaller significances 
in P and r.m.s. scatter for the sSFR split, and the differences 
in r.m.s. are not significant for samples split by Afsteiiar- For 
the a parameters no significant differences are seen. Monte 



Carlo tests on the ^-split sample are similar to those for 

4.3 Effect of assumed cosmology 

Our choice of cosmology could affect our results in two ways. 
The first is the impact on the derived absolute host galaxy 
properties, such as Mstdiar or sSFR. Changing the assumed 
Ho (70kms~^ Mpc~^) has only a superficial effect - all the 
host galaxy masses or other derived properties will change 
relative to each other in the same consistent way. Altering 
the other cosmological parameters, such as JIm, has a more 
subtle effect. In a fiat universe, a smaller f^M will make 
the higher-redshift hosts more massive compared to those 
at lower-redshift (effectively they become more distant); a 
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23 


-0.038 


0.025 






10.90 


56 


-0.037 


0.014 


36 


-0.032 


0.019 


20 


-0.049 


0.021 


Low-z 




7.88 


2 






















9.01 


3 


0.115 


0.031 


















9.79 


6 


0.090 


0.052 


















10.32 


20 


0.018 


0.030 


















10.94 


38 


-0.028 


0.021 














Sample 




Z 




All SNe 






s < 1 






s > 1 






(12- 


hlog(0/H)) 


N 


Mean residual 


Error 


N 


Mean residual 


Error 


N 


Mean residual 


Error 


SNLS 




8.24 


41 


0.023 


0.018 


7 


0.042 


0.035 


34 


0.018 


0.021 






8.56 


35 


0.020 


0.018 


8 


0.070 


0.060 


27 


0.011 


0.018 






8.72 


35 


0.059 


0.019 


8 


0.106 


0.024 


27 


0.035 


0.024 






8.88 


66 


-0.034 


0.017 


35 


-0.020 


0.026 


31 


-0.056 


0.020 






9.00 


15 


-0.022 


0.022 


12 


-0.012 


0.027 


3 


-0.051 


0.033 


Low- 2; 




8.07 


2 






















8.55 


2 






















8.73 


2 






















8.91 


5 


0.049 


0.039 


















9.06 


58 


-0.015 


0.017 















Table 2. Binned SN la residuals for the SNLS and low-redshift samples as a function of host galaxy sSFR, Msteiian a^nd calculated 
with the SiFTO light curve fitter. The host parameter bin centres are given as the mean of all the hosts in that bin. The mean residuals 
are the weighted average of SNe in each bin, and the errors are the errors on that weighted mean. 



larger Q,m will have the opposite effect. Within the errors to 
which cosmological parameters are currently measured, this 
is a small effect - a stellar mass derived with f^M ~ 0.30 
instead of f^M = 0.256 changes by <0.04 dex at z = 0.5 
and <0.07 dex at z = 1. We have run our fits using this 
alternative model and find essentially identical results. 

The cosmological model could also impact the SN prop- 
erties. If the redshift distributions of the host galaxies either 
side of the split point in mass or metallicity are different, 
then any systematic trend in SN brightnesses with redshift 
may increase (or decrease) the significance of the differences 
in the nuisance parameters. For example, an incorrect photo- 
metric zeropoint in one of the SNLS filters would introduce 
a difference in SN la brightness which is a smooth func- 
tion of redshift, and if the mix of hosts also changes with 
redshift, this could introduce a corresponding offset in the 
magnitudes of SNe la in those hosts. 

We test the effect of the assumed cosmology in Fig. [T] 
We start with our default cosmology (^m = 0.256, w = — 1), 
and vary SIm (by ±0.15) and w (by ±0.4) recording the 
difference in derived nuisance parameters. As might be ex- 
pected only small variations in the differences in the de- 



rived nuisance parameters are seen. The largest variation in 
the difference in Mb is with Msteiiar (~ ±0.01), for /3 with 
Mstoiiar (~ ±0.05), and for a with Z (~ ±0.05). Thus the 
assumption of the cosmological model does not drive our 
results. 

4.4 SNLS selection effects 

A final consideration is the possibility of a selection effect 
which operates to either select against fainter SNe la (af- 
ter correction) in massive or low-sSFR galaxies, or brighter 
SNe la in low-Mstoiiar or high-sSFR galaxies. A mechanism 
for the latter is difficult to conceive, but for the former pos- 
sible biases can be imagined. SNe la in high-Afstciiar galaxies 
will be intrinsically fainter (i.e., lower-stretch) and will lie 
against a brighter host background, decreasing their con- 
trast over their host galaxies. If this operates near the SNLS 
spectroscopic limit, then we may be biased to preferentially 
observe the brighter sub-sample of this population in mas- 
sive galaxies, i.e. those that at fixed light-curve shape are 
over-luminous (following the tJint distribution). 

We have partially mitigated against this selection effect 
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Sample 


1( 


3g sSFR 




All SNe 






XI <0 






xi > 










N 


Mean residual 


Error 


N 


Mean residual 


Error 


N 


Mean residual 


Error 


SNLS 




-11.51 


21 


-0.036 


0.029 


15 


-0.045 


0.029 


6 


0.009 


0.093 






-10.03 


37 


-0.021 


0.022 


20 


-0.023 


0.037 


17 


-0.020 


0.021 






-9.31 


61 


0.027 


0.014 


21 


0.042 


0.023 


40 


0.016 


0.018 






-8.39 


70 


0.027 


0.017 


14 


-0.002 


0.040 


56 


0.034 


0.019 


Low-z 




-11.36 


21 


0.017 


0.031 


















-10.06 


31 


-0.007 


0.026 


















-9.38 


20 


0.026 


0.035 


















-8.42 


6 


-0.092 


0.076 














Sample 


loi 


? ■'^stellar 




All SNo 






xi < 






xi > 








(Mq) 


N 


Mean residual 


Error 


N 


Mean residual 


Error 


N 


Mean residual 


Error 


SNLS 




8.31 


25 


0.042 


0.030 


4 


0.041 


0.069 


21 


0.042 


0.035 






9.08 


32 


0.027 


0.019 


8 


0.068 


0.035 


24 


0.012 


0.022 






9.61 


39 


0.063 


0.016 


11 


0.042 


0.029 


28 


0.073 


0.020 






10.27 


32 


-0.017 


0.024 


12 


-0.006 


0.045 


20 


-0.026 


0.026 






10.89 


60 


-0.036 


0.018 


35 


-0.031 


0.024 


25 


-0.045 


0.027 


Low-z 




7.88 


2 






















9.01 


3 


0.074 


0.042 


















9.81 


7 


0.038 


0.063 


















10.30 


21 


0.021 


0.029 


















10.95 


45 


-0.009 


0.023 














Sample 




Z 




All SNe 






xi <0 






xi^O 






(12- 


hlog(0/H)) 


N 


Mean residual 


Error 


N 


Mean residual 


Error 


N 


Mean residual 


Error 


SNLS 




8.25 


38 


0.023 


0.024 


7 


0.043 


0.042 


31 


0.016 


0.029 






8.55 


30 


0.023 


0.019 


8 


-0.047 


0.032 


22 


0.037 


0.022 






8.72 


37 


0.069 


0.017 


10 


0.064 


0.034 


27 


0.071 


0.019 






8.88 


65 


-0.023 


0.018 


32 


-0.014 


0.028 


33 


-0.036 


0.023 






9.00 


16 


-0.037 


0.028 


12 


-0.020 


0.035 


4 


-0.084 


0.040 


Low- 2; 




8.07 


2 






















8.55 


2 






















8.73 


2 






















8.92 


6 


-0.010 


0.062 


















9.07 


66 


-0.000 


0.018 















Table 3. As Table [2] but for the SALT2 light curve fitter. Note that the number of SNe la in each bin can vary from that in Table [2] 
due to slightly different SN parameter culls for SiFTO and SALT2. 



LOG sSFR Low sSFR hosts High sSFR hosts 



split (yr 1) 


NsN 


a 


/3 


Mb 


r.m.s. 


Ns-N 


a 


P 


Mb 


r.m.s. 


-9.30 


96 


1.376±0.128 


2.876±0.159 


-19.193±0.012 


0.127 


100 


1.757±0.181 


3.729±0.162 


-19.126±0.017 


0.143 


-9.50 


74 


1.305±0.148 


2.928±0.184 


-19.196±0.015 


0.123 


122 


1.771±0.161 


3.567±0.148 


-19.130±0.015 


0.151 


-9.70 


61 


1.379±0.148 


2.728±0.231 


-19.208±0.016 


0.120 


135 


1.576±0.149 


3.445±0.136 


-19.149±0.014 


0.147 


-9.90 


53 


1.338±0.170 


2.719±0.232 


-19.204±0.018 


0.123 


143 


1.576±0.145 


3.455±0.136 


-19.152±0.013 


0.145 


-10.10 


27 


1.582±0.255 


2.976±0.451 


-19.205±0.026 


0.128 


169 


1.394±0.114 


3.372±0.125 


-19.167±0.011 


0.142 



Table 4. SiFTO fits with a fixed cosmological model for low and high sSFR galaxies with different values of sSFR (given in column one) 
used to split the sample. The nuisance parameters are free in the fit. 



LOG Afstcllar High-Mstcllar hostS Low-Mstdlar hostS 



split (Mq) 


A^SN 


a 


/3 


Mb 


r.m.s. 


A^SN 


a 


/3 


Mb 


r.m.s. 


9.50 


120 


1.443±0.138 


3.380±0.143 


-19.189±0.013 


0.142 


75 


1.660±0.189 


3.601±0.250 


-19.117±0.020 


0.147 


9.75 


103 


1.501±0.143 


3.258±0.154 


-19.199±0.014 


0.141 


92 


1.612±0.173 


3.644±0.188 


-19.122±0.017 


0.146 


10.00 


92 


1.554±0.148 


3.159±0.163 


-19.206±0.014 


0.141 


103 


1.638±0.164 


3.714±0.170 


-19.121±0.016 


0.143 


10.25 


77 


1.555±0.162 


3.203±0.180 


-19.213±0.016 


0.143 


118 


1.673±0.159 


3.499±0.156 


-19.127±0.015 


0.141 


10.50 


56 


1.577±0.165 


2.937±0.195 


-19.222±0.017 


0.131 


139 


1.515±0.141 


3.480±0.145 


-19.148±0.013 


0.142 



Table 5. As Table g] but for A4tcllar instead of sSFR. 
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Figure 4. As FiglS] but for Afstollar instead of sSFR. 



Z split 






High-Z hosts 










Low-Z hosts 






(12+log(0/H)) 


iVsN 


a 




Mb 


r.m.s. 


Ns-N 


a 




Mb 


r.m.s. 


8.70 


102 


1.512±0.147 


3.238±0.158 


-19.195±0.014 


0.142 


93 


1.661±0.175 


3.800±0.192 


-19.115±0.018 


0.149 


8.75 


92 


1.577±0.152 


3.172±0.165 


-19.204±0.015 


0.143 


103 


1.688±0.161 


3.848±0.176 


-19.110±0.016 


0.146 


8.80 


81 


1.601±0.155 


3.102±0.170 


-19.210±0.015 


0.139 


114 


1.612±0.160 


3.747±0.165 


-19.124±0.015 


0.144 


8.85 


61 


1.534±0.185 


3.095±0.194 


-19.217±0.018 


0.139 


134 


1.739±0.148 


3.614±0.151 


-19.125±0.014 


0.143 


8.90 


34 


1.526±0.187 


2.795±0.225 


-19.218±0.021 


0.113 


161 


1.491±0.130 


3.619±0.139 


-19.151±0.012 


0.148 



Table 6. As Table g] but for Z instead of sSFR. 



by restricting our analysis to those SNe la lying at 2 < 0.85, 
away from the redshift limit of SNLS. The total Malmquist 
bias (including spectroscopi c selection) on our SNe below 
this redshift is < 0.015 mag l|Perrett et al.ll2010l ). compared 
with the size of the magnitude difference in our results of 
^O.lmag. We have also tested for the existence of this ob- 
servational selection effect by examining the SN la residuals 
versus the percentage increase of the SN flux over its host. 
Below a percentage increase of 100% (i.e., the SN is as bright 



as the host measured through a s mall aperture ) , iden tifica- 
tion becomes more difficult (e.g. iHowell et al.l I2OO5I ) . Our 
key diagnostic would be to see brighter SNe (after correc- 
tion) when the percentage increase is <100%, and fainter 
SNe at percentage increases >100%, if this selection effect 
were serious. We show these data in Figure [S] No effect is 
present in our data; only weak trends are present with the 
percentage increase as expected given the weak correlation 

with A/stcllar- 
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8.0 8.5 9.0 0.0 0.2 0.4 0.6 0.8 1.0 

Host galaxy inferred metallicity (12+L0G(0/H)) Relative number 




8.0 8.5 9.0 0.0 0.2 0.4 0.8 0.8 1.0 

Host galaxy inferred metallicity (12+L0G(0/H)) Relative number 

Figure 5. As FiglS] but for metallicity instead of sSFR. 



As a final test, we have also examined the Hubble resid- 
uals of events which were not followed spectroscopically 
by SNLS, but which wer e subsequently loc ated in a offline 
search of the same data l|Bazin et al.ll2009h and are photo- 
metrically similar to SNe la (Bazin et al., in prep.). The 
number of additional small percentage-increase objects at 
z < 0.85 is only 24 (~10% of our total sample), and have 
the same host-dependent trends as those presented in this 
paper. We therefore conclude that this potential selection 
effect cannot drive our results. 



5 DISCUSSION 

5.1 Comparison to other results 

The results of § [4] can be compared to several studies from 
the literature. The result that high-Mstoiiar galaxies host 
the brighte st SNe la after correction was also found by 
iKellv et al.l |2009) at a ~2.5ct significance using a similar 
low-redshift SN la sample to that used here. 



Usi ng galaxy morpholog ies for the same low-redshift 
sample, iHicken et al.l (|2009a ) found weak trends with host 
galaxy type: SNe la in Scd/Sd/Irr galaxies are fainter than 
those in elliptical galaxies after correction {2a difference), 
and have a smaller r.m.s. scatter and lower reddening. As- 
suming that sSFR closely tracks galaxy morphology, we find 
a similar effect in that fainter SNe are located in high sSFR 
galaxies, but we do not reproduce their result that SNe in 
spiral galaxies have a smaller scatter on the Hubble dia- 
gram; in fact our SNLS data suggest the opposite, and at 
low-redshift we see no significant difference in r.m.s. scatter. 
We also see no evidence that SNe in spirals have less redden- 
ing; SNe in spirals appear red der than those in ellipticals. 

Turning to metallicities, Gallagher et al.l (|200S 1 used 
spiral galaxy emission line measurements to directly mea- 
sure log(0/H) for 16 local host galaxies containing SNe for 
which a Hubble residual was available, but found no signif- 
icant trend - their results indicated that brighter SNe after 
correction may l ie in more meta l rich g alaxies, but only at 
90% confidence. iGallagher et al.l (|2008l ) used another small 
sample (17 host galaxies) of local SNe in E/SO galaxies, us- 
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Figure 6. Joint confidence contours (one and two a) in a and /3 (top left), a and Mb (top rigfit), and Mb and /3 (lower left) for fits 
wfiere the SNe are split according to host galaxy sSFR (left panels) and A/gtellar (right panels). Plots with SNe split by Z are similar to 
those for Afstcllar Eind are not shown. The crosses show the best fit values. 
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Figure 7. The effect of the assumed cosmology on the differences in nuisance parameters derived from SNe la located in hosts either side 
of the default split points. The loft-hand figure shows SN la hosts split by A/atoiian the right-hand figure by Z. Either Qm or w is altered, 
with the other parameters held fixed, and the difference in nuisance parameters derived from SNe la in low-Afstciiar a-nd high-Afstdiar 
(or low-Z and high-Z) hosts plotted. 



ing [M/H] spectral measures to show that brighter SNe la 
(again, after correction) were hosted by more metal rich sys- 
tems with ~98% confidence. 



iHowell et all ^200^ ) previously used 55 SNLS SNe fa 
and techniques similar to those in this paper, but found a 
correlation between Hubble residual and inferred host metal- 
Ucity (or Mstdiar) at only ~ 1.3a. These 55 events are also 
included in th is paper, but note that the ffubble residuals of 
iHowell et all (,2009. ') were based on p hotometry, calibrat ion 
and light-curve fitting tools used in lAstier et alj (l2006ll - 
these have all since improved (see discussions in Guv et al] 



2007 : IConlev et al. I l2008l : iRegnault et al] l2009l : IGuv et al. 
20f Of ). The slop e of ffubble residual with inferred log(0/H) 
as measured by IHowell et all (|2009l ') is -O.fO ± 0.07; fitting 
the binned points in the SNLS-metallicity p lot in Fig. B| 
gives a slope of — O.fS ± 0.06, consistent with IHowell et al. 
|2009l ) at ~f .5(7. 

In summary, some similar trends to those presented in 
this paper have been found in the past, but all at quite low 
significance. This is the first dataset where a dependence 
of corrected SN la luminosities on host properties has been 
definitively detected at >3(j confidence. 
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Percentage increase of SN flux over host 

Figure 8. Residuals from the best- fitting cosmology for the SNLS 
SNe as a function of the percentage increase of the events over 
their host galaxies at the time of spectroscopic observation. SNc 
la with s < 1 are shown as open circles, and those with s ^ 1 as 
filled circles. Brighter SNe la after correction have negative resid- 
uals. The red circles are the mean residuals in bins of percentage 
increase, drawn at the mean value of the percentage increase in 
each bin. The error bars on these points represent the bin width 
in percentage increase, and the weighted error on the mean for 
the residual axis. 

5.2 Colour residuals 

The Hubble residuals as a function of colour for SNLS SNe 
la segregated by the sSFR and A/stoiiar of their hosts are 
shown in Fig. |9l with the best-fitting relations over-plotted. 
The (few) objects with red colours could influence the best- 
fitting disproportionately. We repeat the nuisance variable 
fits, each time reducing the maximum allowed colour by 0.05 
and re-examining the fit /3 (Table [7]). Clearly the errors in /3 
increase as the redder SNe are excluded due to the reduced 
baseline in colour. 

The difference between the /3 in low and high sSFR 
galaxies becomes less significant as the reddest allowed SN 
la colour is reduced - most of the difference seems caused 
by the reddest SNe la. Though the colour variation in SNe 
la is not well understood, high sSFR galaxies are likely to 
show a larger range in dust content than low sSFR galaxies, 
and hence display a larger scatter and/or steeper colour- 
luminosity relations for redder SNe la. Removing these red- 
der events is likely to minimise the effect of any difference 
in P between host galaxy types. It may also be the case that 
the colour-luminosity relation of the reddest SNe is domi- 
nated by the effects of dust rather than intrinsic variation 
(in which case (5 ~ Rb), and hence redd e r SNe favour larger 
values of /3. We note that lFolatelli et al] (|2010h actually find 
an decrease in Rv {Rv = Rb — 1) in a sample of low-redshift 
SNe la when including very red SNe la, though these events 
with {B — V) ^ 1 are much redder than any SNLS SN con- 
sidered in this paper. They speculate that this difference in 
Rv in the red dest SNe la could be a result of circumstellar 
dust (see also ' Wang"2005') . 

With smaller maximum allowed colours, the r.m.s. scat- 
ter of the fits of SNe la in each type are also more consistent. 
Under the hypothesis that the reddest SNe la are the most 



affected by dust this might be expected. SN la host galaxies 
are likely to display a ran ge in effective Ry - the Ga lactic 
average is 2.99±0.27 fe.g. iFitzpatrick fc Massal bOOTl ). and 
there is no reason to expect the range present in SN la hosts 
to be smaller. For an E{B — V) of 0.1, this would generate 
an additional scatter of 0.06 mag in SN peak luminosities 
compared to a SN la with a colour excess of zero. 



5.3 SN stretch or host dependence? 

Given the dependence of SN la stretch on host galaxy pa- 
rameters such as Mstoiiar and sSFR (Fig. [2]), the trends of 
SN la brightness with the same parameters (Fig. [3| could 
simply arise from an incomplete stretch-luminosity correc- 
tion, rather than an intrinsic dependence on any third vari- 
able beyond stretch and colour. For example, if the stretch- 
luminosity relation were more complex than a simple linear 
trend but only a linear variable were fit, this could mani- 
fest as a dependence of corrected SN la luminosity on host 
galaxy parameters. 

We investigate this in two ways. We first examine the 
luminosity trends as a function of host galaxy parameters 
segregated by stretch. The SNe la plotted in Fig.[2]are coded 
according to the stretch of the event, and the residuals of 
both s < 1 and s ^ 1 SNLS SNe la are given in Table [2] 
(Table [3] has the SALT2 equivalent). Generally, the same 
trends with Msteiiar, Z and sSFR are seen for high-s and low- 
s SNe la, though clearly the significance is smaller than for 
the entire population given the number of events is reduced. 

We also experiment with a quadratic stretch-luminosity 
term in an attempt to remove the host dependence. We add 
an additional term to eqn. ([T]) of the form 02(5 — l)'^, shown 
in Fig. 1101 The improvement in the quality of the fits with 
this extra stretch term is not significant (Table [8]) and does 
not remove the host-dependent SN la luminosity trends. 
These tests suggest that the trends of SN la luminosity with 
host parameters are largely independent of SN la stretch. 



5.4 Cosmological implications 

Differences in mean SN la properties when split by host 
galaxy properties, which are not removed by corrections cur- 
rently employed in cosmological analyses, could clearly lead 
to systematic errors in cosmological analyses. Observation- 
ally, SNe la in massive galaxies appear brighter than those 
in less massive galaxies - a similar effect is seen when consid- 
ering sSFR, with SNe la in low sSFR galaxies brighter than 
those in high sSFR galaxies. These differences are significant 
at >3cr. 

Evidence for two populations of SNe la with different 
photometric properties is not by itself alarming, as the nui- 
sance variables in any global fit will average to values appro- 
priate for the combination of SNe la. However, any change 
in the mix of SNe la with redshift could introduce a serious 
effect. We measure the change in SN la host Mstdiar or Z 
as a function of redshift (Fig. Ill|) using the SNLS dataset 
unrestricted in redshift range, and the low redshift and R07 
data. We measure the fraction of host galaxies with Mstdiar 
or Z less than the split points defined in earlier sections in 
each redshift bin. We then make a simple linear fit to these 
values as a function of redshift i.e., a + bz, where a and b are 
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Figure 9. Hubble diagram residuals as a function of C for SNLS SNe la before the colour-luminosity relation is applied, split by sSFR 
(left) and A/gteiiar (right), with the cosmological model fixed and the nuisance variables a, beta, and Mg allowed to vary according to 
the type of host. The over-plotted lines show the best fitting relations (a slope of /3). Note that the exact values of the residuals vary 
between the sSFR and Mgteiiar splits as the nuisance variables are allowed to change. As each set of hosts is allowed to have different 
values of Mg, the differences in SN la luminosity between different host types is not present. 



Reddest C Low sSFR hosts High sSFR hosts 



permitted 


A^SN 


a 


13 


Mb 


r.m.s. 


NsN 


a 


/3 


Mb 


r.m.s. 


0.30 


65 


1.377±0.158 


2.730±0.228 


-19.208±0.016 


0.120 


131 


1.573±0.150 


3.444±0.137 


-19.149±0.014 


0.147 


0.25 


65 


1.377±0.158 


2.730±0.228 


-19.208±0.016 


0.120 


130 


1.521±0.150 


3.361±0.142 


-19.152±0.014 


0.145 


0.20 


64 


1.373±0.158 


2.730±0.228 


-19.208±0.016 


0.120 


129 


1.549±0.151 


3.443±0.152 


-19.149±0.014 


0.143 


0.15 


63 


1.435±0.164 


2.947±0.264 


-19.202±0.016 


0.123 


122 


1.614±0.145 


3.214±0.183 


-19.149±0.013 


0.137 


0.10 


61 


1.548±0.170 


3.346±0.328 


-19.190±0.017 


0.129 


113 


1.600±0.148 


3.617±0.257 


-19.136±0.014 


0.140 



Table 7. As Table [J] but testing the effect of the reddest SNe la on the derived nuisance variables. Similar results are found for SALT2 
fits. 



the fit coefficients. To guard against selection effects, we per- 
form tlie fits with and without the low-redshift data where 
the selection effects are quite different to those in the SNLS 
data. In the case of metallicity, the strength of any redshift 
trend is driven by the form of the Mstoiiar-metallicity re- 
lationship assumed - our default relation that evolves with 
redshift leads to a correspondingly larger redshift trend than 
a relation that is assumed constant with redshift (and also 
leads to a decreasing upper metallicity limit for any given 
redshift, as seen in Fig. Ill|l . If a relation with no redshift 
evolution is used, the metallicity-redshift plot becomes es- 
sentially the same as the Matoiiar-redshift plot. 

As expected from a consideration of popular galaxy for- 
mation models, the host galaxies at higher redshift contain, 
on average, less stellar mass and consequently are likely to be 
of lower metallicity. The strength of the evolution is most 
significant with metallicity in the SNLS-fR07 data, but is 
substantially strengthened in Mgteiiar if the low-redshift host 
data are also included. However, including the low-redshift 
data is likely to overestimate the strength of the real physical 
evolution, as the bulk of the low-redshift is strongly biased 
against host galaxies of low stellar mass. It would be ap- 
propriate to include the low-redshift data if the amount of 



evolution likely in current cosmological samples was required 
to be estimated. 

Given the expectation and evidence of evolution, either 
from selection effects or an underlying physical effect, we in- 
vestigate methods for accounting for it in cosmological anal- 
yses, and the effect that it would have if left uncorrected. 

5.4-1 Addttional nuisance variables 

We consider two ways to include additional parameters in 
the fits. The first and most obvious is to include a further 
linear 7x0; term to eqn. ([TJ, where a; is a variable such as 
A^steiiar, Z Or sSFR mcasurcd from the SN la host galax- 
ies. The second is motivated by the idea of two populations 
of SNe la rather than a continuum in properties, and in- 
stead introduces two independent values for Mb in eqn. ([Sjl 
for SNe la located in different types of host galaxy. Explic- 
itly, SNe la located in galaxies with log(Mstciiar) ^ 10.0 
(or Z > 8.8, or log(sSFR) -9.7) are assigned M%, and 
SNe la in galaxies with log(Mstoiiar) < 10.0 (or Z < 8.8, 
or log(sSFR) > —9.7) are assigned Mb- In principle, dif- 
ferent values of a and /3 could also be fit for SNe la lying 
either side of the split point. However, the results of § 14.21 
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Figure 11. The change with redshift in SN la host galaxy Afstcllar (left) and inferred Z (right). The mean values in bins of redshift are 
shown as red circles, drawn at the mean redshift in each bin, and with error bars representing the bin width for the redshift axis, and 
the error on the mean in the Mstj-iiar a-nd ^ axis. The horizontal dashed lines show the split points used in §|4] Open circles refer to the 
low redshift or R07 samples, and filled circles the SNLS sample. Note that the differences between the low-redshift and SNLS hosts is 
unlikely to be a real physical effect, and is generated by selection effects in the low-redshift sample. The lower panels show simple linear 
fits to the fraction of SN la hosts with log(Mi3toiiar)lOA^0 (or Z < 8.8), including (dashed line) and excluding (solid line) the low-redshift 
data (shown as an open circle). The fit parameters are shown; the numbers in italics include the low-redshift data. 



E m 
E 




Figure 10. Hubble diagram residuals as a function of s for SNe la 
split by Mstoiiar before the stretch-luminosity relation is applied. 
The over-plotted lines show the best fitting linear (solid) and 
quadratic (dashed) stretch-luminosity relations (Table [Sj. Only 
a single Mb is used, and so the SNe in high-A/gtciiar and low- 
-A/stcllar hosts are slightly offset. 



show that these nuisance variables are more consistent than 
Mb, i.e. the Mb captures most of the host galaxy dependent 
variation. 

The results of these different fits with the three different 
host variables are given in Table [8] (again, the cosmology is 
held fixed). We also hold (Tint fixed, so the addition of an ex- 
tra parameter into these fits leads to a significant reduction 
in x^i particularly for the Mstoiiar and related metallicity 
fits. _F-tests show that the addition of an extra parameter 
improves the quality of the fits beyond that expected by the 
addition of the extra parameter. For the two Mb fits, the 
significance is 3.9-4.4(7 for Mstciiar and 4.1-4.4(7 for Z. The 
linear fit, though still improving the x^y generally performs 
more poorly: 3.3-3.5(7 for Afstciiar and 2.6-3.0(7 for Z. Fits 
for two Mb using Afstciiar or Z perform better than those 
with sSFR; note that the linear fits in sSFR are difficult to 
implement as a zero SFR cannot be represented in log(sSFR) 
and must be assigned a fixed nominal value. 

In the SNLS plus low-redshift fits with two Mb, M% 
is more poorly measured than Mb. This can be traced to 
the biased nature of the low-redshift sample. In SN la only 
fits, the absolute magnitudes are heavily influenced by the 
low-redshift data where the SN la brightness is less depen- 
dent on the cosmology, and as there are far fewer SNe la in 
the low-redshift sample in the low-Afstciiar/low-Z group, the 
statistical precision of M% is correspondingly reduced com- 
pared to MJj where there are abundant low-redshift SNe la. 



© 0000 RAS, MNRAS 000.[Hf23l 



SN la host galaxies 19 



Sample 



7 or 02 



Mb or Ml, 



Ml 



DOF 



Base fits: 
SNLS 

SNLS+low-2 
SNLS, two a 
SNLS+low-5;, two a 



1.384±0.103 
1.353±0.083 
1.349±0.107 
1.346±0.083 



3.348±0.124 
3.127±0.100 
3.347±0.124 
3.121±0.101 



1.094±0.943 
0.572±0.733 



-19.170± 0.010 
-19.178± 0.008 
-19.160± 0.011 
-19.172± 0.008 



190.8 
270.4 
189.8 
270.0 



192 
261 
191 
260 



A4tcllar fits: 

SNLS, 7 
SNLS, two Mb 
SNLS+low-2, 7 
SNLS+low-2, two Mb 



1.553±0.115 
1.571±0.115 
1.515±0.093 
1.508±0.091 



3.365±0.127 
3.382±0.127 
3.150±0.103 
3.178±0.103 



0.0423± 0.0125 
0.0412± 0.0106 



-19.172± 0.010 
-19.209± 0.015 
-19.176± 0.009 
-19.211± 0.011 



-19.127± 0.015 
-19.132± 0.014 



179.9 
176.0 
256.4 
251.4 



191 
191 
260 
260 



Metallieity fits: 
SNLS, 7 
SNLS, two Mb 
SNLS+low-2, 7 
SNLS+low-2, two Mb 



1.521+0.116 
1.587+0.116 
1.482+0.092 
1.520+0.092 



3.398+0.129 
3.377+0.127 
3.184+0.104 
3.172+0.103 



0.1124+ 0.0412 
0.1107+ 0.0328 



-19.180+ 0.011 
-19.214+ 0.015 
-19.181+ 0.009 
-19.212+ 0.011 



-19.127+ 0.015 
-19.132+ 0.014 



184.3 
175.2 
261.0 
250.9 



191 
191 
260 
260 



sSFR fits: 
SNLS, 7 
SNLS, two Mb 
SNLS+low-2, 7 
SNLS+low-2, two Mb 



1.580+0.115 
1.525+0.113 
1.541+0.094 
1.470+0.090 



3.254+0.127 
3.295+0.125 
3.078+0.102 
3.107+0.101 



-0.0438+ 0.0118 
-0.0434+ 0.0101 



-19.176+ 0.010 
-19.210+ 0.017 
-19.181+ 0.009 
-19.211+ 0.013 



-19.144+ 0.013 
-19.152+ 0.012 



178.4 
182.2 
253.7 
259.6 



191 
191 
260 
260 



Table 8. Fits involving extra nuisance parameters based on host galaxy measurements, specifically M^^^n^,., Z and sSFR. In all cases 
the cosmology is held fixed. 



The use of two absolute magnitudes has other advan- 
tages over a hnear host parameter term. There is no theoret- 
ical reaso n to suppose that a ny dependence on metallieity 
be linear: pTimmes et al.l (|2003l l show that the putative effect 
will be strongly non-linear. From a practical viewpoint, the 
use of two absolute magnitudes is substantially less sensitive 
to systematics in host galaxy parameter determination. The 
categorisation of low-Mgteiiar versus high-A/gtciiar is not sen- 
sitive to systematic errors associated with the choice of base 
SED library or assumptions about IMF or dust extinction, 
nor to the statistical errors on the masses, particularly at 
low-Msteiiar where these can be considerable. 

The other nuisance parameters are quite stable to the 
introduction of a 7 or M% term. (3 is almost unchanged, and 
Q tends to drift to larger values by about 0.15. This latter 
effect can be understood when considering the stretch de- 
pendence on host galaxy properties. Though the trends in 
SN la brightness are the same for both low and high stretch 
SNe 15. 3p . lower-s SNe la are more prevalent in massive 
and low sSFR galaxies. These SNe la prefer a brighter AIb 
than those in low-Afstoiiar hosts, making their Hubble resid- 
uals, when fitting for two Mb, more positive (fainter). This 
in turn will lead to a systematic steepening of the stretch- 
luminosity relation, which is reflected in the larger values of 
a seen when using third parameters. 

5.4.2 Systematics in cosmology 

What is the size of the effect of this host galaxy dependence 
in cosmological analyses? We analyse this by fitting the low- 
redshift, SNLS and R07 data (314 SNe la) with two simple 
cosmological models with and without the host galaxy term, 
and assessing the variation of the fit parameters. For sim- 
plicity, we use only the statistical errors of the SNe - in a full 
cosmological analyses, the covariance matrix between differ- 



ent SNe accounting for all systematic errors should be used 
IConlev et al.ll201oi ): here we only aim to assess the size of 
the effect. Our cosmological models are a flat, ACDM model 
(fitting only for f^M, i.e. assuming a constant w = —1), and 
one in which w is still constant but not required to be equal 
to —1 (fitting for Q,m and w). In each case the crint for each 
SN sample is adjusted to give a reduced of 1.0. We also 
show the effect of adding external constraints into the fit 
using the Baryon Acoustic Oscillation (BAO) constraints of 
lEisenstein et al.l (|2005l ). The fit results, expressed as a shift 
in the cosmological parameters from their values without 
the host parameter, can be found in Table [5] We use the 
cosfitter progran^j, written by one of us (AC) and ad- 
justed to fit for two Mb, to perform the fits. We test both 
A/stoiiar and gas-phase Z as the third variable on which to 
split the data. 

The size of the effect - the shift in the cosmological pa- 
rameters with a third variable - is comparable to the simple 
statistical precision (Table l^}. These shifts are larger when 
metallieity is the third variable (> la) instead of Msteiiar 
(< la). One source of systematic error is the choice of split 
point in the host variable. We investigate the dependence of 
the cosmological results on this choice in Fig. 1121 We show 
the effect on w (including BAO constraints) and the as 
the split point is varied. The median error in mass is ~ 0.14 
dex (~ 0.03 once converted to gas phase Z), and the relative 
change in w as the split points are varied over this range is 
small. One approach would be to iterate on the split point in 
cosmological fits in order to find the model with the smallest 

As a final point, we note that the difference in /3 be- 
tween low-sSFR and high-sSFR host galaxies, and a change 



^ http : //qold. astro .utoronto . ca/conley/simple_cosf itter/ 
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Third ACDM w model, SN only w model with BAO 

variable SUm (ct) SQm (f^) Si" (rr) 5Qm (o") (ct) M]j M| 

A^teiiar 0.012 (0.6) 0.068 (0.8) 0.173 (0.8) -19.12±0.03 -19.19±0.02 0.001 (0.1) 0.040 (0.6) -19.13±0.02 -19.20±0.02 

Z 0.031 (1.6) 0.066 (0.8) 0.205 (1.0) -19.10±0.03 -19.19±0.02 0.006 (0.3) 0.078 (1.2) -19.11±0.02 -19.20±0.02 



Table 9. Simple effect of a host galaxy parameter on cosmological constraints. The shift in the cosmological parameters are shown when 
a host parame ter is added, togethe r wit h the shift expressed as a fraction of the simple statistical error-bar. A full cosmological analysis 
is presented in lConlev et al. 1 1 I2OIOI) and lSuUivan et al. I 1I2OIOI ). 




5.5 Comparison to theory 



10 

LOG M,.„, (Ms) 




i.4 8.6 

Z3„„[12+log(0/H)l 



9.0 9.2 



Figure 12. The sensitivity of the cosmological results on the 
choice of split point used for the host galaxy variable. Cosmologi- 
cal fits are performed at each choice of split point, and the fraction 
of SNe la with the host parameter less than the split point (upper 
panels), the of the fit (middle panels), and the change in w 
(lower panels) are plotted as a function of the split point. The up- 
per plots show Af stellar I the lower ones Z. The default values arc 
shown as vertical lines. The w measures include BAO constraints. 
When all SNe la are assigned to one of the groups as happens on 
the left and right edges of these plots, the cosmological results are 
identical, crint is held fixed in these tests at the values required 
to produce a reduced of one without additional parameters - 
the decrease in can be interpreted as a decrease in (Tint when 
including the host parameter. 



in the mix of the two with redshift, will introduce a mild 
redshift evolution in the mean value o f /3. Such an evolution 
has been claimed (jKessler et al.ll2009l ). tho ugh is not presen t 
to any significant extent in SNLS3 data (|Guv et al.ll2O10l ). 
Given the small evolution in Mstciiar or sSFR with redshift in 
the SNLS data, we estimate any /3 evolution would be <0.1 
across our entire redshift range, below the detectability of 
our current dataset. 



iTimmes et al] l|2003h showed that the amount of ^®Ni gen- 
erated in a SN la explosion will depend on the metallicity of 
the progenitor star. More metal-rich stars will generate in- 
creased ^^N during hydrogen burning as C, N and O are con- 
verted to in the CNO cycle. During core helium burning, 
this is converted to neutron-rich ^^Ne, which will favour 
the production of stable, neutron-rich ^*Ni at the expense 
of the radioactive ^^Ni that drives the luminosity of SNe la. 
Fainter SNe la should therefore be produced in more metal- 
rich environments. 

While at first sight this may appear to disagree with 
the results in Figs . [4] and [5] which show the opposite trend, 
iKasen et al.l (|2009l ) show this metallicity efTect will also af- 
fect the SN light-curve duration - SNe la in metal-rich envi- 
ronments are both fainter and have narrower light curves (or 
smaller stretches). Though this trend operates in the same 
sense as the global stretch-luminosity relation, and is there- 
fore partially corrected by it, the correspondence is not exact 
and a difTerence is expected - for a fixed light-curve width, a 
higher metallicity SN la will be brighter. Though the slope 
of the metallicity-driven stretch-luminosity relation is ex- 
pected to change with metallicity, over the stretch range that 
our SNLS data sample (80% of our SNe have 0.9 < s < 1.2 
or 1.2 > Amis > 0.8 usin g the relation from lConlev et al] 
(l2008t )'). lKasen et all (|2009l ) show that the metallicity should 
act with a fairly constant gradient in stretch-luminosity 
space. The effect of metallicity should therefore translate 
to a simple offset in the stretch-luminosity relation derived 
from our data, rather than a change in slope. This is equiva- 
lent to a change in Mb, with brighter SNe in metal-rich en- 
vironments. Thus this theoretical expectation of progenitor 
metallicity on SN la luminosity appears qualitatively con- 
sistent with our data. (Note that these models also predict 
that SNe la in massive galaxies should have slightly smaller 
stretches due to their high metallicities, all other variables 
being equal.) 

Clearly, other physical variables correlate with host 
properties such as Mstdiar and sS FR, including stel lar age 



L g Stel la 

and dust content (see discussion in lKellv et al.l l200!l). Dust 
appears less likely to be responsible for the effects in the 
SNLS data. The colour-luminosity relations of SNe la in 
low- Mstciiar aud high-Msteiiar hosts are consistent, as are the 
mean colours of the SNe in those hosts, as well as the distri- 
bution of those colours. There is also the expectation that 
SNe la will be affected by only small amounts of extinction, 
particularly in flux lim ited surveys like SNLS (|Hatano et al.l 
1 19981 : ICommindliooi ) . 

The role of the progenitor system stellar-age is sub- 
stantially more complex, and clearly cannot be ruled out 
in our data - the most massive hosts will on average con- 
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tain stars of older ages. At the time of writing, no clear 
prediction has been made as to the effect of progenitor age 
on SN la luminosity. There is speculation that progenitor 
age drives the global relationship between light-curve shape 
and host galaxy type, where morphologically earlier-type 
galaxies host lower-stretch SNe la. Of importance in this 
context is the result that the Afstdiar-SN-brightness depen- 
dency does not depend on the stretch of the SN. Low-s 
SNe in lower- Mgteiiar hosts have a similar corrected bright- 
ness as high-s SNe in low-A/stciiar hosts, both consistently 
fainter than their counterparts in high-Matoiiar hosts. Under 
the assumption that age is the major driver of the stretch- 
luminosity relationship, this would suggest that it is metal- 
licity that is responsible for the SN la Mstoiiar-luminosity 
dependence that we observe. 

If metallicity is the key physical variable driving the lu- 
minosity trends, then the existence of radial me tallicity gra- 
dients in galaxies fe.g. llfenrv fc Worthevl[l999l ) will compli- 
cate the analysis. Studies generally show that the metallicity 
of galactic disks decreases with increasing galactocentric ra- 
dius, implying that a simple segregation of SNe based on 
host galaxy stellar mass or metallicity may (ultimately) re- 
quire refinement to include, for example, the location of the 
SN within the host. No trends are apparent between SN 
residual and galactocentric radius in the SNLS data. How- 
ever, many factors complicate this analysis, the most perni- 
cious of which is that the measure d SN positions must be 
deprojected fe.g. Ilvanov et a l."2000''). This requires a knowl- 
edge of the structural parameters of the galaxy (for example 
inclination), difficult to measure for high-redshift galaxies 
using ground-based data. Additionally, for SNe la with long 
delay times between progenitor formation and explosion, the 
progenitor star may migrate away from the region in which 
it is formed. Investigations of this type must await higher- 
quality low-redshift data. 



6 CONCLUSIONS 

In this paper we have examined the photometric properties 
of Type la Supernovae (SNe la) as a function of their host 
galaxy properties using new data from the Supernova Legacy 
Survey (SNLS) and literature data for SNe la at low-redshift 
and z > 1. Our principal findings are 

• As expected, SN la light-curve widths closely track 
the specific star-formation rate (sSFR) and stellar mass 
(Mgteiiar) of their host galaxies - massive and/or low sSFR 
galaxies host SNe la with lower stretches (narrower light 
curves). There is only a mild dependence of SN la colour 
(C) on the host galaxy - SNe la in low sSFR galaxies are 
slightly bluer than those in high-sSFR host galaxies, but the 
colours of SNe la in high-Mstoiiar and low-Msteiiar hosts are 
consistent. 

• SNe la in low sSFR host galaxies, and SNe la in massive 
host galaxies, are systematically brighter by 0.06-0.09 mag 
(> 3a") after light-curve shape and colour corrections. This 
is not dependent on any assumed cosmological model. In- 
terpreting Mstoiiar as an approximate metallicity indicator, 
this implies that metal-rich environments host brighter SNe 
la post-correction. 

• Galax;ies with low sSFRs host SNe la with a smaller 
slope (/3) between SN la luminosity and colour, and which 



have a smaller scatter on Hubble diagrams. This is depen- 
dent on the maximum colour of the SN which is considered 
- the effect is more significant (2-3a") when redder SNe are 
included. The slope of the relationship between stretch and 
luminosity is consistent for SNe la across different galaxy 
types. 

• The SN la luminosity variation with host properties in- 
troduces a systematic error into cosmological analyses, as 
the mean Mstoiiar and sSFR of the hosts, and hence the mix 
of SNe la, evolves with redshift. This effect is exacerbated 
by the biased selection of existing low-redshift SNe la, and 
amounts to a systematic error on w comparable to the sta- 
tistical errors. 

• We propose that this be corrected for in all future SN la 
cosmological analyses by incorporating a host galaxy term 
into the fit. Specifically, we show that the use of two absolute 
magnitudes for SNe la, one for those in low-A/stciiar (or low- 
metallicity) hosts, and one for events in more massive, metal- 
rich hosts, leads to an improvement in cosmological fits at 
3.8-4.5cr, and removes the host dependence. 

• The SN la luminosity effects appear consistent with 
theoretical expectations of the dependence of SN la lumi- 
nosities on progenitor metallicity. The effect of metallicity is 
predicted to translate into an offset in the effective stretch- 
luminosity relation, with brighter SNe in metal-rich environ- 
ments after correction, as observed in our data. 

These results are consistent with earlier observational 
studies which have found evidence for SN la luminosity 
variation as a function of host galaxy properties at 2- 
2.5a confidence ([Gallagher et all l2008l : iHicken et all l2009bl : 
iKellv et al.l \200^ ). The SNLS dataset. spanning a larger 
range in host Afatciiar and a less biased selection of SNe la 
than at low redshift, provides tighter constraints on host 
galaxy dependencies and broadly supports the results of 
these earlier studies. 

Redshift evolution in SN la properties driven by 
metallicity effects in the SN la population has been sug- 
gested as a systematic error in S N la cosmology for many 
vears jHoflich et al.l Il998l I2OO0I: IPommguez et all l200li : 



iTimmes et al.ll2003l : ]Kasen et al.l l2009''). The results of this 
paper detect this evolution, demonstrate that it is qualita- 
tively consistent with theoretical predictions, and show that 
it can be corrected using supplementary information on the 
environments in which the SNe explode. There are two key 
implications from the perspective of measuring dark energy. 
The first is that additional information on the SN la host 
galaxies, such as multi-colour photometry covering a broad 
wavelength range, will be an essential requirement for fu- 
ture SN la cosmological analyses and surveys. Multi-colour, 
rolling searches similar to SNLS are obviously well-placed to 
provide this information as a natural product of the survey 
strategy. 

The second implication is the urgent need for new low- 
redshift SN la samples, where events are selected with- 
out regard to their host galaxy type or brightness. Exist- 
ing searches are mostly galaxy targeted - repeatedly imag- 
ing the same catalogued galaxies - and result in heav- 
ily biased sa mples of SNe la. The Nearby Supernova Fac- 
tory ( Aldering et al.l l2002t ) . and the next generation low- 
redshif t surveys such a s the Palomar Transient Factory 
(RTF; iRau et all l2009l : iLaw et al.l |2009| ) and SkyMapper 
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l|Keller et al.ll2007l ). should provide samples of SNe la se- 
lected in a similar way to those at higher redshift, reducing 
the systematic uncertainties associated with evolving galaxy 
populations. 
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